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Abstract. We summarise recent progress in the understanding of the rotational evolution of 
low-mass stars (here defined as solar mass down to the hydrogen burning limit) both in terms 
of observations and modelling. Wide-field imaging surveys on moderate-size telescopes can now 
efficiently derive rotation periods for hundreds to thousands of open cluster members, provid- 
ing unprecedented sample sizes which are ripe for exploration. We summarise the available 
measurements, and provide simple phenomenological and model-based interpretations of the 
presently-available data, while highlighting regions of parameter space where more observations 
are required, particularly at the lowest masses and ages > 500 Myr. 
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1. Rotation period data 

A compilation of most of the available rotation period (and some v sini) measurements 
in open clusters for stars with masses M < 1.2 M Q is shown in Figure Q] (a list of 
references is included in Table [1]). We plot rotation period as a function of stellar mass, 
rather than using the more conventional quantities of colour or spectral type on the 
horizontal axis, since the diagram spans ages from the early pre-main-sequence (PMS), 
at ~ 1 Myr in the ONC, to the main sequence (MS), and neither colour nor spectral 
type are invariant for a given star over this age range. The masses are, of course, model- 
dependent, but the majority of binning in mass and morphological examination used 
in this work require only that the mass scale is approximately correct, which should be 
reasonably well-satisfied by the PMS stellar evolution tracks. We use those of the Lyon 
group, from Baraffe et al. (1998) and I. Baraffe (private communication), throughout. 

An expanded version of Figure [T] omitting many of the more sparsely sampled clusters 
and all the vsini observations, is shown in Figure [21 By examining the evolution of the 
morphologies of these diagrams, we can already draw some (model-independent) conclu- 
sions regarding the evolution of the rotation periods in these clusters. For simplicity, we 
discuss two broad mass ranges: stars close to solar mass (e.g. 0.9 < M/M Q < 1.1), and 
fully-convective stars (M < 0.4 M Q ), which represent the "tail" of the distribution which 
emerges especially in the NGC 2547 and NGC 2516/M35 panels of Figure[2] 

Considering the solar mass stars first, it is clear that the behaviour on the PMS (which 
lasts to ~ 30 Myr for a solar mass star, until it reaches the zero- age main sequence, here- 
after ZAMS) and on the MS up to the age of the Hyades, varies as a function of rotation 
rate. The slowest rotators in the ONC have periods of ~ 10 days. The most basic predic- 
tion we can make for the evolution of the rotation rate is to assume angular momentum 
is conserved, in which case as the star contracts towards the ZAMS, it should spin up. 
However, for these slowly-rotating stars, the period remains approximately constant all 
the way to NGC 2362, and has spun up to - 8 days by the age of NGC 2547 (~ 40 Myr), 
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Figure 1. Compilation of 3100 rotation periods (and some wsini measurements) for stars with 
masses M < 1.2 Mq in young (< 1 Gyr) open clusters, from the literature. Plotted in each panel 
is rotation period as a function of stellar mass for a single cluster. The appropriate references for 
each panel are given in Table Q] All the masses used in t his contribution were computed using 
the I-band luminosities of the sources and the models of Baraffe et al. (1998) | assuming values 
of the age, distance modulus and reddening for the clusters taken from the literature. 



at which point the contraction ceases. Over the same age range, the most rapid rotators 
gradually evolve from a period of ~ 1 day at the ONC, to ~ 0.6 days at NGC 2362 
(~ 5 Myr), and ~ 0.2 days at NGC 2547. This follows closely the prediction from stellar 
contraction. We therefore conclude that there is some mechanism removing angular mo- 
mentum from only the slow rotators, to an age of ~ 5 Myr, and a short time after this 
the angular momentum losses cease, leaving the star free to spin up for a short time until 
it reaches the ZAMS. The net result of this is to yield a wide range of rotation rates in 
the early-MS clusters such as the Pleiades and M35. 

By the age of the Hyades however, the rotation rates converge onto a single well- 
defined sequence. It is clear that this process is well underway by the age of M37 (~ 
485 Myr). There must therefore be another rotation-rate-dependent angular momentum 
loss mechanism operating on the early-MS, such that more rapid rotators lose more 
angular momentum, to drive all the stars toward the same rotation rate at a given mass. 

In contrast, the behaviour of the lowest mass stars is clearly rather different. On the 
PMS, these stars all appear to spin up rapidly, and reach very rapid rotation rates at 
the ZAMS (for a 0.4 M Q star, the ZAMS is reached at - 150 Myr), and furthermore, 
the maximum rotation period seen is a very strong function of mass (this is most clear 
in the NGC 2516 diagram in Figure [2]). This indicates that there is little rotation rate 
dependence in this mass domain due to the similarities between the morphologies of the 
diagrams especially from NGC 2264 to NGC 2516, and furthermore, that these stars 
lose little angular momentum on the PMS, in contrast to the slowly rotating solar mass 
stars. Moreover, between NGC 2516 and Praesepe, the limited quantity of data available 
suggest that there is essentially no evolution of the rotation period, so little angular 
momentum appears to be lost on the early-MS, again in contrast to the solar-type stars. 
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Figure 2. Enlarged version of Figure[T] showing a few selected clusters which have large samples 
of rotation periods available. NGC 2516 and M35 have been combined into a single panel, since 
they have essentially identical ages, as have the Hyades and Praesepe. 



We proceed in fj2]to examine the physical interpretation of these observational results. 



2. Models of rotational evolution 

Recall from |JT] that we needed to invoke two mechanisms of angular momentum re- 
moval: one operating for ^5 — 10 Myr on the PMS, and one operating on the main 
sequence, with the property that it must produce a convergence in the rotation rates 
between ~ 100 and 600 Myr. 

2.1. Pre-main- sequence angular momentum losses 

For the first of these, the PMS angular momentum loss, the ~ 5 — 10 Myr timescale 
suggests an obvious candidate: the circumstellar discs that surround these stars in the 
earliest stages of their evolution dissipate on these timescales (e.g. Haisch et al. 2001). 
The precise mechanism by which angular momentum is removed due to the presence 
of a disc is still a matter of debate, with the currently-favoured hypothesis being an 
accretion-driven stellar wind (e.g. Matt & Pudritz 2005), or "disc locking" (e.g. Konigl 
1991; Collier Cameron et al. 1995). For our purposes, we shall simply treat the effect 
of the disc as holding the angular velocity of the star constant (by removing angular 
momentum) for the lifetime of the disc. 

In reality, there will be a range of disc lifetimes, and thus the stars will be "released" 
from this angular velocity regulation at a range of times, giving rise to a spread of rotation 
rates on the PMS and ZAMS, as required by the observations. Moreover, examining one 
of the early-PMS clusters shows us a "snapshot" of this process in action. The best- 
observed example is the ONC, shown in Figure [H 

If we presuppose that stars begin with a tight distribution of initial rotation periods, 
centred around 8 — 10 days, in a histogram of rotation period we should expect to see a 
population of stars at this period, and then a tail of faster-rotating stars, corresponding 
to the ones now uncoupled from their discs, with faster rotating stars having uncou- 




Figure 3. Distribution of rotation periods for ONC stars with masses M > 0.25 Mq (using the 
D'Antona & Mazzitelli 1998 models; this corresp onds to M > 0.4 M for the Baraffe et al. 1998 



models used in this contribution), replotted from Herbst et al. (2001) The two panels show the 
same distribution binned in linear period (left) and log period (right). 

pled earlier. This is indeed what is seen in Figure [3J where it is instructive to plot the 
distribution in log period rather than the more conventional linear period. 

Moreover, this hypothesis makes an observationally-testable prediction: the population 
of slowly-rotating stars should show measurable indications of the presence of a disc (e.g. 
mid-IR excess) or of active accretion (i.e. be classical T-Tauri stars, CTTS), whereas the 
rapidly-rotating stars should not have discs, or have recently-dissipated discs, and not 
be active accretors (i.e. weak- lined T-Tauri stars, WTTS). 

In reality, extracting this observational signal from the cluster data has proved to be 



difficult (see Rebull et al. 2004 for a detailed discussion of this issue). Rebull et al. (2006) 
presented the first statistically significant detection of this effect, using mid-IR excess 
measurements from the IRAC instrument aboard the Spitzer space telescope, the current 
method-of-choice, finding that slowly-rotating stars arc indeed more likely to posses discs 
than rapidly-rotating stars. However, in doing so this study revealed a wrinkle in the 
argument: they also found a puzzling population of slowly-rotating stars without discs. 
Several subsequent studies (e.g. Cieza & Baliber 2007) have confirmed these results, 
which now appear to be placed on a firm footing. 

2.2. Main sequence angular momentum losses 

The mechanism for angular momentum loss on the main sequence, at least for solar 
mass stars, is thought to be a solar-type magnetised stellar wind. The time-dependence 
of rotation rates in this age range has been firmly established observationally for solar- 



type stars since Kraft (1967) and from the age of the Hyades to the Sun, obeys well 
the famous Skumanich (1972) t 1 / 2 law (i.e. u> oc t~ x l 2 ). This can be reproduced on a 
more theoretically motivated framework from parametrised angular momentum loss laws 
(usually based on Kawaler 1988; the t 1 ^ 2 law can be reproduced by setting n = 3/2 and 
a = 1 in his model). 



Although the Skumanich (1972) law provides a good description of the evolution of 
solar-type stars from the age of the Hyades to the Sun, and for the majority of stars 
in the Pleiades to the age of the Hyades, simply evolving the Pleiades distribution for- 
ward in time to the age of the Hyades using this law produces stars which are rotat- 
ing much more rapidly than observed in the Hyades. Therefore, most modellers mod- 



ify the Kawaler (1988) formalism to incorporate saturation of the angular momentum 
losses above a critical angular velocity u sat (Stauffer & Hartmann 1987; Barnes & Sofia 
1996). The saturation is further assumed to be mass-dependent, to account for the mass- 
dependent spin-down timescales observed on the early-MS. 
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The resulting angular momentum loss law assumed in our models is: 

dt -\ „ ( R\ 1/2 ( M^ 1 ' 2 {2A) 



We note that although saturation is reasonably well-motivated observationally (e.g. by 
the saturation observed in the relationship between X-ray activity and rotation rate), this 
is a rather unsatisfactory feature of the present models, since there is a certain amount 
of arbitrariness in the way this parameter is introduced (for example, the choice of power 
of ui in Eq. 12.11 for the saturated angular velocity dependence of dJ/dt is arbitrary). An 
important area for future work is to develop a theoretical framework for this phenomenon, 
and for solar- type winds in general. 

2.3. Comparison of models to data 



The models of rotational evolution we use are described in detail in Bouvier et al. (1997) 



features in this contribution. 



Allain (1998) Irwin et al. (2007) and Irwin (2007) We summarise only the most salient 



The methodology of Bouvier et al. (1997) which we adopt here essentially implements 



a rotational evolution model by using a standard non-rotating PMS stellar evolution 
track to compute the variation of the stellar parameters as a function of time, principally 
the moments of inertia of the radiative and convective regions of the star, which are then 
used as an input to determine the angular velocity as a function of time. 

This model in the simple form described has four parameters (which could all be 
functions of mass): the initial angular velocity u>q, normalisation of the solar- type angular 
momentum loss law K (see Eq. 12. ip . saturation velocity w sa t, and the lifetime of the 
circumstellar disc Tdi sc - For solar- type stars, we can fix K by requiring that the model 
reproduces the observed rotation rate of the Sun. 

Figure [4] shows the result of fitting this model to the rotation period data using a 
simple nonlinear least squares routine. For the moment, we have required the radiative 
core and convective envelope to have the same angular velocity, i.e. a "solid body" model. 

The results show that the model does a reasonable job for the late-time evolution 
from the Hyades to the Sun. but this is largely by construction since we used the 



Skumanich (1972)] law as an input! At earlier times, in general, we find that the rapid ro- 
tators are better-reproduced by the model, with the only major issue being at ~ 150 Myr 
where the model underpredicts the upper bound to the observed rotation rates. For the 
slow rotators, the model rotates too rapidly given the ONC and NGC 2362 as an initial 
condition around the ZAMS and early-MS. The disc lifetimes behave as expected, being 
shorter for the rapid rotators. 

The difficulty fitting the slow rotators motivates the introduction of an additional 
parameter into the model. In particular, we now relax the assumption of solid body 
rotation, and allow the core and envelope to have different angular velocities, coupling 
angular momentum between them on a timescale r c (this is done in a fashion which 
tries to equalise their angular velocities). r c = represents the solid body case already 
considered. The net effect of this "core-envelope decoupling" from the point of view of the 
observed rotation rate (which is that of the envelope) is to "hide" angular momentum in 
the core when it forms, which is then coupled back into the envelope gradually, providing 
a "late-time replenishment" of the surface angular velocity. 

Figure [5] shows the result of applying this revised model to the solar-type stars. We can 
see that the fit is substantially improved. Importantly, we find that the value of r c for the 
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Figure 4. Rotational angular velocity u) plotted as a function of time for stars with masses 
0.9 < M/Mq ^ 1.1. Crosses show the open cluster rotation period data such that each cluster 
collapses into a vertical stripe on the diagram, and short horizontal lines show the 25th and 
90th percentiles of u, used to characterise the slow and fast rotators respectively. The lines show 
rotational evolution models for 1.0 Mq stars, fit to the percentiles for each cluster using a simple 
unweighted least squares method. For this plot, we have assumed the stars rotate as solid bodies 
(i.e. constant u as a function of radius inside the star). Plotted are the ONC, NGC 2264, NGC 
2362, IC 2391, IC 2602, a Per, M35, M34, the Hyades, and the Sun. 



slow rotators is relatively large, ~ 110 Myr indicating inefficient core-envelope coupling, 
whereas for the fast rotators it is short, ~ 6 Myr: these stars have efficient core-envelope 
coupling and rotate more like solid bodies. A corollary of this is that slow rotators develop 
a higher degree of rotational shear across the convective/radiative boundary. 



This result may have important observationally-testable consequences. Bouvier (2008) 



discusses one of these: the impact on the surface abundance of lithium. In particular, 
the higher rotational shear in the slow rotators is likely to induce additional mixing in 
these stars, which may bring Li down into the star where it can be burnt. The net effect 
would be a higher level of lithium depiction in slow rotators than fast rotators, and 
these should therefore display a lower lithium abundance on the main sequence. Indeed, 



Soderblom et al. (1993)| report that this appears to be the case in the Pleiades. 

Having examined the models for solar-type stars, we now proceed to the fully-convective, 
very low- mass stars (M < 0.4 M©). Figure [6] shows the results of applying the same mod- 
els to this mass domain. Since these stars are fully-convective, there is no possibility for 
core-envelope decoupling, so they should rotate as solid bodies. 

The model is clearly a poor match to the observations. In particular, the solar-type 
wind appears to loose too much angular momentum at late-times, resulting in a sharp 
drop in the predicted rotation rates starting at ~ 200 Myr, which does not appear 
to be seen, although the available data are still very limited at present, with vsini 
measurements giving only upper limits for the slowest rotating stars. This cannot be 
"fixed" by invoking saturation, since this has already been done for the model shown in 
Figure [6l which has w sat = 1.8 ujq (i.e. every star is saturated in the relevant age range). 
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Figure 5. As Figure [4] but now allowing the radiative core and convective envelope of the stars 
to have different angular velocities, coupling angular momentum between the two on a timescale 
t c . We have also omitted NGC 2264, since the models presently have difficulty explaining the 
observations in this cluster. 

It would not be surprising that a solar-type wind does not reproduce the evolution of 
fully-convective stars, since this is thought to be launched at a tachocline of shear at the 
radiative/convective interface, and these stars posses none. However, as highlighted by J. 
Stauffer in the questions, we should be cautious in overinterpreting the rather sparse data 
presently available at these masses and ages, especially given the large number of wsini 
upper limits and the natural bias of the rotation period method toward more active, 
and hence presumably more rapidly-rotating, stars, resulting from the extremely small 
photometric amplitudes of the rotational modulations at these old ages. 

This domain of very low-masses and late (> 500 Myr) ages represents an important 
area for future observational studies to explore, and could provide us with a valuable 
insight into the winds of fully-convective stars. 
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